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ABSTRACT 

Mergers of double neutron stars are considered the most likely progenitors for short gamma-ray bursts. 
Indeed such a merger can produce a black hole with a transient accreting torus of nuclear matter (Lee & 
Ramirez- Ruiz 2007, Oechslin & Janka 2006), and the conversion of a fraction of the torus mass-energy 
to radiation can power a gamma-ray burst (Nakar 2006). Using available binary pulsar observations 
supported by our extensive evolutionary calculations of double neutron star formation, we demonstrate 
that the fraction of mergers that can form a black hole - torus system depends very sensitively on the 
(largely unknown) maximum neutron star mass. We show that the available observations and models 
put a very stringent constraint on this maximum mass under the assumption that a black hole formation 
is required to produce a short gamma-ray burst in a double neutron star merger. Specifically, we find 
that the maximum neutron star mass must be within 2 — 2.5M Q . Moreover, a single unambiguous 
measurement of a neutron star mass above 2.5Mq would exclude a black hole - torus central engine 
model of short gamma-ray bursts in double neutron star mergers. Such an observation would also 
indicate that if in fact short gamma-ray bursts are connected to neutron star mergers, the gamma-ray 
burst engine is best explained by the lesser known model invoking a highly magnetized massive neutron 
star (e.g., Usov 1992; Kluzniak & Ruderman 1998; Dai et al. 2006; Metzger, Quataert & Thompson 
2007). 
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for massive stars (Belczynski et al 2008). The code em- 
ploys state-of-the-art predictions for neutron star and 
black hole masses based on hydrodynamic core collapse 
simulations (Fryer & Kalogera 2001) and detailed stel- 
lar structure and evolution calculations for massive stars 
(Timmes, Woosley, & Weaver 1996). Our models predict a 
Galactic NS-NS merger rate in the range ~ 10 — lOOMyr -1 
(Belczynski et al. 2007a), in good agreement with the em- 
pirical estimate of ~ 3 — 190Myr _ (Kim, Kalogera, & 
Lorimer, 2006). The spread in our predicted rates orig- 
inates from including the most significant model uncer- 
tainties associated with the treatment of dynamical mass 
transfer episodes (common envelope phases) , which are in- 
volved in the formation of most double compact objects 
(Belczynski et al. 2007a). 

2. RESULTS 

In Figure Q] we compare short GRB rates with NS-NS 
merger rates in the present-day (redshift 0) universe. Ex- 
trapolating the NS-NS merger rates to the local universe 
by assuming a star-forming density of 10 -2 Milky Way- 
equivalents per Mpc 3 El, we estimate the local universe NS- 
NS merger rate to be in the range ~ 100 — 1000Gpc~ 3 yr _1 . 
By comparison, the estimated conservative lower limit on 
the short GRB rate is ~ 10Gpc~ 3 yr -1 , based on the 
BATSE/SWIFT sample (Nakar 2006). This estimate re- 
lies on very conservative assumptions: (i) there is no colli- 
mation and (ii) there are no bursts dimmer than we have 
already observed, thus providing a true lower limit on the 

1 Introduction of different modes of star formation and redistribution of progenitors over old elliptical and younger spiral galaxies does not 
increase the predicted merger rates by more than a factor of ~ 3 O'Shaughnessy et al. 2007. 
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1. INTRODUCTION 

Gamma-ray bursts (GRBs) have been separated into 
two classes: long-soft bursts, and short bursts (Nakar 
2006,Gehrels et al. 2007). The origin of long-soft bursts 
has been connected to the death of low-metallicity massive 
stars (Piran 2005, Gehrels et al. 2007). However, while ob- 
servations support a binary merger origin for short bursts 
(Nakar 2006, Gehrels et al. 2007), the exact nature of the 
progenitor remains uncertain: they could be either dou- 
ble neutron stars (NS-NS) or black hole - neutron star 
(BH-NS) binaries. The number of BH-NS binaries that 
both merge and produce GRBs is hard to estimate since 
(i) no such system has yet been observed; (ii) formation 
models are rather uncertain and predict very small BH- 
NS merger rates (likely too small to explain most of the 
short bursts); and (hi) theory suggests that the fraction of 
BH-NS mergers producing bursts depends sensitively on 
the black hole spin and spin-orbit orientation (Belczyn- 
ski et al. 2007b), but black hole birth spins are not well 
constrained observationally or theoretically. On the other 
hand, NS-NS binaries are only observed in the Milky Way, 
but their properties and numbers are also in agreement 
with theoretical models, and their merger rate is sufficient 
to explain the present-day short burst population (Nakar 
2006, Belczynski et al. 2007a). 

We have performed an extensive theoretical study of 
high-mass binary stars (potential progenitors of NS-NS 
systems) using StarTrack, a population synthesis code in- 
corporating the most up-to-date and detailed input physics 
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rate. Therefore, even adopting the most optimistic predic- 
tions for the NS-NS merger rate and the most pessimistic 
bound on the local short GRB event rate, the fraction of 
NS-NS mergers f„ r b that produce GRBs must be greater 
than at least 10 to explain the majority of known short 
bursts. 

In this paper we start with the assumption that all short 
GRBs are connected with NS-NS system mergers that pro- 
duce a black hole. We discuss the implications of relaxing 
this stringent assumption at the end of the paper. 

From our models we also derive physical properties of 
double neutron stars, with individual masses of neutron 
stars being of particular interest. Figure [2] shows the rela- 
tion between progenitor (single star) mass and final rem- 
nant mass used in our evolutionary calculations. Mass 
transfer and other binary interactions change this simple 
picture, through both accretion and mass loss, which can 
either increase or decrease an individual binary compo- 
nent mass. However, (Belczynski et al. 2007b) argues 
that we do not expect significant mass accretion onto the 
components of NS-NS binaries. The population model 
we adopt for our discussion here produces NS mass distri- 
butions that appear consistent with the current observed 
NS-NS sample, at least in the extent of the mass ranges 
(Fig. 3). While mass transfer does influence the remnant 
masses (e.g., smearing the narrowly peaked mass distribu- 
tion implied by Fig. 2), the qualitative structure is largely 
preserved, as one would expect from isolated stellar evo- 
lution combined with an initial mass function that falls 
steeply with increasing initial mass. The predicted neu- 
tron star mass distribution only very weakly depends on 
evolutionary model assumptions because the neutron star 
formation mass is almost independent of the progenitor 
mass (Timmes et al. 1996) and mass accretion in NS- 
NS progenitor binaries is rather small (Belczynski et al. 
2007a). 

Depending on the masses in the progenitor binary and 
the highly uncertain nuclear equation of state, the final 
remnant of a NS-NS merger may or may not collapse to 
a black hole. We estimate the final gravitational mass of 
the compact remnant as 

M rem = 0.9(M nM + Af ns , 2 - 0.1M Q ), (1) 

where the initial neutron star masses are denoted by 
M ns .i, M nSi 2 and we have assumed that the torus mass 
is sufficiently large to power a GRB (i.e., ~ 0.1M©) (Seti- 
awan, Ruffert, & Janka 2004, Lee & Ramirez-Ruiz 2007) 
and that 10% of rest mass is lost in neutrinos (Lattimer 
& Yahli 1989; Timmes et al. 1996). Higher rest mass 
loss (e.g., Metzger, Thompson & Quataert 2007) would 
only strengthen our subsequent conclusions. Because stars 
more massive than 18M (progenitors of massive neutron 
stars with M ns ~ 1.8M Q ) are much rarer than those form- 
ing lighter neutron stars (M ns ~ 1.35M Q ) we a priori ex- 
pect that most remnants from NS-NS mergers will have 
rather low mass M rem ~ 2.3M Q (see eq. 1 for two 1.35M Q 
neutron stars). 

Neither observations nor nuclear theory have yet pinned 
down the maximum neutron star mass Afjs^aj above 
which a black hole must form. Thus, the fraction of binary 
neutron stars which produce black holes and are able to 
power short GRBs is set by the fraction of mergers such 
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We therefore calculate the fraction of our simulated NS- 
NS mergers that lead to black hole formation and a short 
GRB as a function of M nSjlnax ; see Fig. 4. Observations 
of the highest mass neutron stars (<; 2M Q ; ( Barziv et al. 
2001, Ransom et al. 2005) and lowest mass black holes (;> 
3M Q ; (Orosz 2003, Casares 2006) only weakly constrain 
this parameter. Remnant masses from NS-NS mergers 
(AT rcm ) obtained both from our simulations and from ob- 
servations all fall very close to the range 2.2 — 2.5M©. 

Comparing Figs 1 and 4 we immediately deduce that, 
since the fraction / gr b of NS-NS mergers that produce 
short GRBs must be greater than 10 -2 (Fig. 1), the neu- 
tron star maximum mass M ns , max must be less than 2.5M Q 
(Fig. 4) . Because we lack a robust lower bound on the mass 
of the residual torus surrounding the black hole, we have 
adopted a conservative upper limit on Af nS:max obtained 
by assuming a negligible torus mass (i.e., replace 0.1M Q 
with in eq. 1). 

This result has been obtained with the assumption that 
all (k = 1.0) short GRBs are connected with NS-NS merg- 
ers. It is however possible that only a fraction of short 
GRBs is produced in NS-NS mergers. How does our result 
depend on the fraction k of short GRBs that are connected 
with NS-NS mergers? The lower limit on / gr b is then de- 
creased by k, see Figure [T] If k J> 0.1, then the limit lower 

limit becomes / gr b > 10~ 3 , and as is clearly seen from Fig- 
ure [3l the upper limit on the maximum mass of a neutron 
star remains unchanged. For the values of k <J 0.01 the 
NS-NS mergers are not important for overall short GRB 
population, as the mergers would consist of only 1% of 
the short GRBs. Therefore, in this case short GRBs do 
not provide information about the merger product. 

3. DISCUSSION 

Our proposed limit on the maximum neutron star mass 
is still above the maximum masses allowed by almost all 
proposed models for the nuclear equation of state (Lat- 
timer & Prakash 2007). However, the proposed limit 
would remain unchanged even if a dramatic improvement 
in short GRB surveys led to a significantly larger lower 
bound on the local short GRB rate: because of the sharp 
decrease in / gr b with M nSjmax shown in Fig. 4. If, however, 
electromagnetic observations could constrain the least lu- 
minous short GRBs and thus provide an upper bound 
on the short GRB rate, with gravitational wave obser- 
vations at the same time accurately determining the NS- 
NS merger rate, then / gr b could also be constrained from 
above. If only a fraction of NS-NS mergers produce short 
bursts, because / gr b depends so sensitively on M nSjmax , the 
combination of upper and lower limits would constrain the 
maximum neutron star mass extremely tightly, even if the 
assumptions going into eq. 1 are relaxed. 

While our limit on the effective maximum neutron star 
mass is entirely empirical, detailed merger models in- 
cluding realistic relativistic dynamics, neutrino transport, 
magnetic fields, and potentially even energy extraction 
from the final black hole remain under intense investi- 
gation (Janka & Ruffert 1996; Oechslin & Janka 2006). 
Many merger remnants are expected to be (temporarily) 
rotationally supported against collapse (Morrison, Baum- 
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garte & Shapiro 2004), with a "hypermassive" remnant 
neutron star eventually spinning down and collapsing to a 
black hole (Faber et al. 2006; Duez et al. 2007; Shibata & 
Taniguchi 2006). Our model only relies upon the current 
consensus on double neutron star mergers, as summarized 
by Oechslin & Janka (2006): sufficiently massive binary 
mergers produce a black hole and only mergers that pro- 
duce a black hole extract enough energy to power short 
GRBs. 

Known Galactic black holes extend in mass up to 
10 — 15M Q (Casares 2006), while two recently discovered 
black hole candidates in other galaxies (Orosz et al. 2007, 
Prestwich et al. 2007) have even higher masses of ~ 16 
and ;> 2AM q. Clearly black holes can form with rather 
high masses in different types of environments. The lower 
mass limit is not well constrained observationally, as the 
highest-mass neutron stars barely reach 2M Q , while the 
lowest- mass black holes are above 3M . In order to ex- 
plain the observed short GRBs with NS-NS mergers, un- 
der the assumption of black hole - torus central engine 
model, we have shown that the maximum neutron star 
mass must be lower than 2.5M Q . However, pulsar surveys 
(Ransom et al. 2005) have discovered increasingly more 
massive neutron stars. So far in our analysis we have in- 
cluded only NS-NS systems formed in the field. There is 
one known relativistic double neutron star system in the 
Galactic globular cluster M15, that has probably formed 
through dynamical interactions. This binary consists of 
two low-mass neutron stars (1.36 and 1.35Mq; Jacoby et 
al. 2006) very similar to those in the Galactic field, so the 
results of our analysis are not changed by this isolated ob- 



servation. Moreover, it was estimated that no more than 
10 — 30% short GRBs can originate from mergers of dou- 
ble neutron stars formed in globular clusters (Grindlay, 
Portegies Zwart, & McMillan 2006). 

If any observation can be made that establishes unam- 
biguously a pulsar mass (either in the field or a globular 
cluster) over 2.5M Q , this would exclude a black hole - 
torus short GRB central engine model for double neutron 
star mergers. We note that a tentative mass measurement 
for a pulsar of 2.74 ± 0.21M Q was recently reported by 
(Freire et al. 2007). If this measurement is confirmed, the 
double neutron star mergers may still be possible progen- 
itors for short GRBs. However, the central engine model 
will need to be reexamined. In particular, it was proposed 
that a merger of two neutron stars may lead to the forma- 
tion of a magnetar; a rapidly rotating highly magnetized 
and high mass neutron star (with or without a torus) that 
can lead to a short GRB (e.g., Usov 1992; Kluzniak & Ru- 
derman 1998; Dai et al. 2006; Metzger et al. 2007). Grav- 
itational wave observatories (LIGO, VIRGO) may provide 
the direct evidence that NS-NS merger can produce a short 
GRB if there is a coincidence of the burst and the inspiral 
gravitational wave signal. It may also be possible to dis- 
tinguish a merger product (NS versus BH) from the shape 
of the merger and ringdown signal or from radio pulses if a 
magnetar was formed in a nearby gamma-ray burst event. 

We would like to thank Neil Gehrels, Duncan Lorimer, 
Scott Ransom, Ben Owen, Jerome Orosz, Chris Stanek, 
John Beacom, Paulo Freire, Chunglee Kim, Todd Thomp- 
son and an anonymous referee for very useful discussions. 
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Fig. 1. — Comparison of the double neutron star merger rates and short GRB event rates. The solid black line and 
arrows indicate a firm lower bound on the short GRB event rate (Nakar 2006), based solely on the rate of detected bursts. 
Depending on the amount of beaming and the fraction of distant faint short GRBs that are missed, the true event rate 
is often estimated to be at least 10 times larger (Nakar 2006). This lower limit is smaller than the double neutron star 
merger rate estimated for the Milky Way both from (i) observations of Galactic binary pulsars (filled blue region) and 
(ii) our population synthesis simulations (filled red region), when these two estimates are extrapolated to cosmological 
scales. Based on the maximum plausible double neutron star merger rate with the minimum plausible short GRB event 
rate, the fraction / gr b of binary mergers that lead to short GRBs should be greater than 10 -2 if double neutron stars are 
the progenitors of short GRBs. 
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Fig. 2. — Initial (Zero Age Main Sequence) mass to final compact object mass relation for single stars. This represents 
our current understanding of compact object formation. Stars below about 7.5M form white dwarfs; stars in the narrow 
range around 8M can potentially form very light neutron stars through electron capture supernovae (Podsiadlowski et 
al. 2004). More massive stars show a well defined bifurcation caused by different modes of energy transport in the stellar 
core: stars below 18M form light neutron stars (~ 1.35M ), while stars above this mass form heavy neutron stars 
(~ 1.8M Q ). Above ~ 2OAf stars experience partial fallback of material that can turn nascent neutron stars into black 
holes. Compact objects originating from stars of ~ 20 — 22Af form either very heavy neutron stars or low-mass black 
holes depending on the unknown limiting mass between these two remnant types (expected to lie around 2 — 3M ). 
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Fig. 3. — Predicted mass distribution for neutron stars in merging double neutron star binaries. First born neutron 
stars are slightly heavier as they can accrete some matter from their unevolved binary companions. Population synthesis 
models (red and blue lines) are shown along with measured neutron star masses for the known double neutron star 
binaries. Although more observations are needed to constrain the shape of this distribution, the mass ranges of observed 
and predicted systems are in agreement. We use direct mass estimates for B1913+16, B1534+12, J0737— 3039 and 
J1756— 2251 (O'Shaughnessy et al. 2008), while for J1906+0746 we assume that both neutron stars have masses of 1.3Mq 
(total system mass is 2.6AT ; Lorimer et al. 2006). The few compact objects found in our simulations with masses as 
high as ~ 2.5M© may well be low- mass black holes (see also Fig. [TJ. 
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Fig. 4. — Gamma Ray Burst production efficiency as a function of the maximum neutron star mass in the framework 
of the double neutron star model for GRBs involving the formation of a black hole. The mass of the merger product is 
plotted here as the blue (observations) and red (theory) lines. There is a sharp drop in number of NS-NS systems that 
can form a merger with mass over 2.5M@. For example, only 1 in 10 3 NS-NS mergers can form a remnant with a mass of 
2.5M Q or higher. Therefore, if short GRBs are connected to NS-NS mergers, the maximum neutron star mass is required 
to be M ns , max < 2.5Mq. For comparison we show observed masses of the lowest-mass black holes (GRO J0422— 32, GRS 
1009—45; Casares 2006) and highest-mass neutron stars (Vela X-1, Tcrzan 51, and Terzan 5J; Barziv et al. 2001, Ransom 
et al. 2005). These observations, along with our findings, constrain the maximum mass of a neutron star to lie in the 
narrow range of 2 — 2.5M Q . 



